Glass Transition for Driven Granular Fluids 
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We investigate the dynamics of a driven system of dissipative hard spheres in the framework of 
mode-coupling theory. The dissipation is modeled by normal restitution, and driving is applied 
to individual particles in the bulk. In such a system, a glass transition is predicted for a finite 
transition density. For increasing inelasticity, the transition shifts to higher densities. Despite the 
strong driving at high dissipation, the transition persists up to the limit of totally inelastic normal 
restitution. 



PACS numbers: 64.70. P-, 64.70.ps, 61.20.Lc, 64.70.Q- 

The celebrated jamming diagram of Liu and Nagel [l[ 
generated a lot of interest in recent years. The conjecture 
is that in the space spanned by the parameters packing 
fraction tp, temperature T, and external stress a, there 
is a region where the material is solid like or jammed, cf. 
Fig. [T] This perspective unifies the concepts of jamming 
of macroscopic, athermal particles and of the glass tran- 
sition of microscopic, thermal particles. A lot of work has 
been devoted to the point J — the arrest of static or quasi 
static granular assemblies Q. Similarly, the glass transi- 
tion of a supercooled molecular liquid has been studied 
extensively, corresponding to a transition line in the tp- 
T-plane, which in the case of hard spheres is parallel to 
the T-axis. In such a system of elastic hard spheres, the 
glass transition is described well by the mode-coupling 
theory (MCT) [Hi] . Whether or not the unified picture 
of the jamming diagram holds is still a matter of debate, 
though. 




FIG. 1: a) The jamming diagram for sheared systems as a 
function of inverse density tp~ , temperature T, and shear 
stress a. b) Jamming diagram for driven inelastic hard 
spheres where the driving power v^ T replaces the shear stress. 
In this latter case, the temperature dependence is trivial for 
T > 0, and the origin of the graph lies at random-close pack- 
ing (rep). 

It was shown by theory and computer simulation that 



both Newtonian and Brownian equations of motion yield 
the same glassy dynamics 0-0]. The situation is differ- 
ent when the system is subject to shear. MCT has been 
extended recently to the case of colloidal suspensions un- 
der shear Q, and for this case it was shown that any 
finite shear rate is able to destroy the glass transition. 
While remnants of the glass transition are still affecting 
the dynamics, full arrest is no longer possible. Another 
scenario is proposed by the above mentioned jamming 
diagram where applied shear stress unjams the system 
but can be compensated by higher density or lower tem- 
perature, cf. Fig.QTa). 

It is the objective of the present paper to investigate 
the possibility of a glass transition in a granular fluid in 
a steady state when dissipation is balanced by bulk driv- 
ing. Experiments by Abate and Durian [9| and Reis et al. 
[1 01 ] showed indications of a granular glass transition in 
such fluidized granular systems in two dimensions. These 
observations would fit into a modified jamming diagram, 
where the shear-stress axis a is replaced by an axis that 
quantifies the driving force t>Dr used to compensate for 
the dissipative interactions among the granular particles. 
In the following, we shall investigate if a granular glass 
transition can exist and how such a transition can be de- 
scribed by an appropriate theory. It will be demonstrated 
that (1) the combination of MCT Q with granular kinetic 
theory predicts a glass transition for a driven dissi- 
pative system, (2) the nature of the transition depends 
on the degree of dissipation, and (3) granular dynamics 
cannot be scaled onto Brownian or Newtonian dynamics. 

We consider the non-equilibrium stationary state of 
a driven granular fluid comprised of N — > oo identical 
hard spheres of diameter d and mass to in a volume V. 
Energy dissipation in binary collisions is modeled by a 
coefficient of normal restitution e. Due to the energy 
loss in the collisions the system needs to be driven in 
order to achieve a stationary state. We use a simple bulk 
driving mechanism, e.g., as in air fluidized beds [9j. The 
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particles are kicked with frequency WDr according to 

v' i {t)=v i (t) + vn I £ i (t). 

The driving amplitude, t>Dr, is constant and the direc- 
tion of the kick, is chosen randomly from a Gaus- 
sian distribution, P(£), with unit variance. To ensure 
momentum conservation, we choose pairs of neighboring 
particles and kick them in opposite directions [ljj]- A 
stationary state is reached, when the energy loss due to 
collisions is balanced by the energy input due to driving. 
A simple estimate relates the properties of the driving 
(uiDi.fDr) to the collision frequency u> co n- 

T 

W co u(l -£ 2 )- = WDr^Dr • 

The time evolution of the system consists of ballistic 
motion in between binary collisions and random kicks. 
These can be formally incorporated in a Pseudo-Liouville 
operator C+ [III EH i which generates the time evolution 
of an observable, such as the density 

p i^> = ~jy 5Z exp (* 9 ' ri (*)) = cxp(itC + )p q (Q) . 

i 

For the discussion of the long-time dynamics, the central 
quantity of interest is the density correlation function, 

F(q,t) = J d 3 6P(0 J dTw(T)p* q (0)p q (t) =: ( Pq (0)\p q (t)) , 

which is directly accessible from computer simulations 
and experiments. Here, w(T) is the stationary iV-particle 
distribution and V = {r^ Vi}f =1 denotes a point in phase 
space. 

Following Mori and Zwanzig, the normalized correla- 
tion function, <j> q (t) = F(q,t)/S(q) = F(q,t)/F(q,0), can 
be represented in terms of restoring forces and a memory 
kernel 14|, 



(d 2 + u q d t + n 2 q ) 4> q {t) 



dTM q (t-T)d T <t> q (T\ (1) 



with 

N N 2 



The thermal velocity vq = \JT jm and initial conditions 
are <f> q (0) = 1,9*^(0) = 0. 

Detailed balance does not hold and consequently, the 
transition rates of forward and backward reactions are 
not simply related, (p q \£ + j q ) ^ (j q \£, + p q ) , as it would 
be the case in equilibrium systems. For the same reason, 
the memory kernel, M q (t) — (i?J|i? 9 (t)), is now given by 
the cross-correlation of two unequal fluctuating forces, 
R q 7^ R q , driven by the reduced dynamics. Details of 
the calculation can be found elsewhere The repre- 
sentation in Eq. ([T]) is exact. It correctly accounts for 



the conservation laws for the density, p q (t) and the lon- 
gitudinal momentum, qjh{t) — 9tp q (t). Energy is not 
conserved in a granular medium and transverse momen- 
tum is decoupled. Hence the representation guarantees 
the correct hydrodynamic limit of <\> q (t) . 

To proceed, we have to resort to approximations. First, 
we need an approximate form of the iV-particle distri- 
bution to compute static (equal time) correlations. We 
assume that positions and velocities are uncorrelated 
w(T) = w r ({ri})w v ({vi}) and that the velocity distri- 
bution factorizes into a product of one-particle distri- 
butions, w v ({vi}) = Y\.i w li v i)- The precise from of 
w\{v) is not needed, it only has to satisfy (v) = and 
(v 2 ') = 3T/m < oo. Furthermore, the system is assumed 
to be isotropic and homogeneous except for the excluded 
volume: w r ({r.i}) — Yli<j^i r ij ~ For the restoring 
forces we find 
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J (qd) + 2j 2 {qd)} , 
e 1 — e 



with uje the Enskog frequency for the elastic case. 

An additional approximation is necessary to compute 
the memory kernel. The success of MCT for the descrip- 
tion of dense molecular and colloidal fluids motivates its 
application also to dissipative granular fluids. First, we 
project the fluctuating forces onto products of densities. 
Second, the resulting higher order correlations are fac- 
torized into pair correlations. Thereby, one generates an 
explicit expression for M q (t) = m q (t)fl 2 in terms of <fi q (t) 



m q [4>](t) « Me)—r I Sk V qk <f>k(t)<f> q -k(t) (2] 



with V q k given by 

V q k = SkS q -k [q ■ kc k + q ■ (q - k) c q - k f 

The direct correlation function, c q , is simply related to 
the static structure factor via the Ornstein-Zernike equa- 
tion nc q = and A q {e) = [1 + (1 -e)S q /(l + e)] _1 
depends on e explicitly. Inserting the mode coupling 
approximation, Eq. ©, into Eq. ([T]), we get a self- 
consistency equation for the scattering function <fi q (t). 
The only further input that is required is the static struc- 
ture factor S q . For simplicity, we use the elastic Percus- 



Yevick expression [14[ here. Future work [15( will study 
the influence of a more precise structure factor that de- 
pends on the coefficient of restitution e. 

A memory function under driving is not necess arily 
positive and might not even be a real function [16j . 
Hence, it is surprising that for a driven granular fluid 
the only change compared to the elastic case is the e- 
dependent prefactor A q (e) > 0, see Eq. Conse- 
quently, the memory kernel itself is positive and a vast 
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amount of work that has been devoted to the mathemat- 
ical structure of the standard mode coupling equations 
(for a compilation see, e.g., p) is readily applicable to 
the granular system. This finding is quite remarkable 
because it implies that a large number of results derived 
for equilibrium systems also holds for a system far from 
equilibrium. In particular, a positive memory function 
in Eq. guarantees positive spectra. On the other 
hand, the well established universality of glassy dynam- 
ics is broken for granular fluids since the memory kernel 
depends on e explicitly. 
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FIG. 2: Transition density ip c as a function of the coefficient 
of restitution e. The short horizontal bar indicates the result 
for the elastic case for e = 1. The inset shows the evolution 
of the critical exponent a with e. 

A glass transition is signaled by the appearance of a 
time persistent part of the density correlations f q :— 
limt-j.oo (j) q {t) in Eq. (Q]). In this limit, Eq. ([T]) reduces 
to the algebraic equation f q /(l — f q ) = m q [f] which is 
solved readily by standard procedures [l7j]. For all val- 
ues of the coefficient of restitution, < e < I, an ideal 
glass transition of the driven granular fluid is found with 
a transition density <p c (e), cf. Fig. [3J For increasing 
dissipation, i.e., smaller e, the glass transition is shifted 
to higher densities. This can be understood as follows: 
For increased dissipation, e < 1, the prefactor A q {e) in 
Eq. ([2]) becomes smaller than unity. Dissipation and driv- 
ing hence weaken the memory effects and destabilize the 
glass. This needs to be compensated by a higher density. 
The resulting compensation of enhanced dissipation and 
driving by increased density can be represented in a cor- 
responding jamming diagram as shown in Fig. [TJ), where 
for simplicity the mean field relation oc (I — e) 2 /4 
was used. It might be reassuring that, although the crit- 
ical density increases with increasing dissipation, it stays 
below the density of random-close-packing. 

The full dynamics of Eq. ([1]) is shown in Fig. [3J The 
time dependent scattering function 4> q (t) with a density 
close to the glass transition density (f c shows the generic 
two step relaxation. After an initial fast relaxation, the 
scattering function approaches a plateau 4> q (t) — f q and 
decays to zero only for very long times, provided the den- 
sity is still below the critical one. Close to the transi- 




FIG. 3: Dynamics of the coherent density correlator g (t) for 
the wave vector qd = 4.2 at the respective critical densities 
V? c (e) and at ip — 0.999ip c (e) for e = 1.0 (dashed lines, elastic 
case), 0.5 (full lines), and 0.0 (dotted lines). 

tion point a critical decay is found onto these plateaus, 
4> q {t) — f q oc t~ a . The variation of the critical exponent 
a is shown in the inset of Fig. [2] For values below the 
transition, ip < p c , a second power law describes the de- 
cay from the plateau, known as the von-Schweidler law, 
<f) q (t) — f q oc t b . The exponent b (not shown here) also 
varies with the coefficient of restitution e, and is uniquely 
related to a. 




qd 

FIG. 4: Critical glass form-factors /| for coefficient of restitu- 
tion e — 1.0 (empty squares, elastic case), 0.5 (filled circles), 
and 0.0 (filled diamonds) as a function of wave number qd. 

At the transition point and beyond, ip > tp c , the cor- 
relation function assumes a finite long-time limit f q >0 
which sets in at a critical plateau value These values 
are shown in Fig. 0] for various values of e. The increasing 
dissipation has three noticeable effects: (1) correlations 
at small wavenumbers are enhanced, (2) oscillations re- 
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fleeting the local structure become less pronounced, and 
(3) the localization length (indicated by the inverse of the 
width of the f q distribution) decreases. The last finding 
is a consequence of the glass transition taking place at a 
higher density, cf. Fig. [2J 

It was shown earlier that Newtonian (N) and Brownian 
(B) systems show the same glassy dynamics 0- This is 
shown for the MCT dynamics in Fig. [5] for curves N and 
B, where N only needed to be shifted along the time axis 
to match B. In contrast, due to the explicit dependence 
on e in Eq. (|2]h for any s < 1, the granular long-time 
dynamics (G) cannot be scaled on top of the Newtonian 
or Brownian results. In addition, for granular dynamics 
at different e there also exists no single master curve, cf. 
Fig. El Therefore, granular dynamics leads to a funda- 
mentally different long-time behavior. 
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FIG. 5: Dynamic scattering function 4> q (t) for qd = 4.2 at 
the glass transition, ip — ip c , and for slightly lower volume 
fraction, ip = 0.999<£ c for Newtonian dynamics (N) with v q = 
0, Brownian dynamics (B), and granular dynamics (G) with 
coefficient of restitution e = 0.5. The Newtonian dynamics is 
scaled along the time axis to match the Brownian dynamics 
at long times. No such rescaling is possible for the granular 
dynamics. The inset shows three lines of equal diffusivity 
(from top to bottom: D = 0.7, 0.8, and 0.9 in relative units) 
taken from the data of a computer simulation fl2] ]. 

In conclusion, we have shown that MCT can be ex- 
tended to granular fluids, which are in a steady state 
far from equilibrium. It can be shown that the result- 
ing memory kernel in Eq. (J5|) is positive, and this al- 
lows that most mathematical theorems from equilibrium 
MCT carry over to the driven granular case. An ideal 
glass transition is observed for all values of the dissi- 
pation or equivalently all values of the driving. Since 
the balance between dissipation and driving implies that 
the driving amplitude is proportional to the inelastic- 



ity Up r cx (1 — e) 2 , the glass transition defines a line in 
the 95 _1 -versus-u 2 )r plane of a generalized jamming dia- 
gram, cf. Fig. [T] The universality known for Newtonian 
and Brownian dynamics is broken in the granular case 
- however, the predicted differences in the critical ex- 
ponents, cf. inset of Fig. [3J and glass form factors, cf. 
Fig. HI are relatively small. Comparably large changes 
with increased dissipation are expected in the transition 
densities as shown in Fig. [21 This prediction can be sup- 
ported by looking at a precursor of the glass-transition 
line in computer-simulation data: For increasing dissipa- 
tion, points of equal diffusivity are found to be shifted to 
higher densities (inset of Fig. [SJ in accordance with our 
predictions for the glass transition in Fig. [5J Therefore, 
we expect our results to be testable in further computer 
simulation studies and also in experiments. 
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